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Partial-wave and helicity operators for the scattering of two hadrons in lattice QCD 
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Department of Physics and Maryland Center for Fundamental Physics, 

University of Maryland, College Park, MD 20742, USA 

Partial-wave operators for lattice QCD are developed in order to facilitate the identification of the 
spins of two-hadron scattering states corresponding to zero total momentum. Taking the periodic 
boundary conditions for lattice states into account, orthogonal sets of partial-wave operators for 
orbital angular momentum are identified. When combined with the intrinsic spins of the hadrons, 
orthogonal sets of parent operators for total angular momentum J and projection M are obtained. 
The parent operators are subduced to irreducible representations of the octahedral group in order 
to obtain descendant operators for use in lattice calculations. The descendant operators retain 
orthogonality with respect to J. The spin of a state can be identified by the spin of parent operators 
that dominate creation of the state. For nonzero total momentum, operators are developed for a 
range of helicities and they are subduced to irreducible representations corresponding to the different 
directions of total momentum. Sets of operators that include a sufficient range of helicities allow 
identification of spin J when a state couples to operators with helicities less than or equal to J, but 
not to operators with higher helicities. 


I. INTRODUCTION 

Lattice operators subduced from irreducible represen¬ 
tations (irreps) of group, SU( 2), to irreps of the octahe¬ 
dral group, Oh, have been used to determine the spins 
as well as the energies of single-hadron states in lattice 
QCD. Recent works used such operators at zero mo¬ 
mentum to determine the excited-state spectroscopy of 
mesons, based on two-quark operators |T], and baryons, 
based on three-quark operators [1]. 

Subduction from 517(2) associates each lattice opera¬ 
tor with values of J and M , where J is the total angular 
momentum and M is its projection along the z-axis. 

Using operators that are smeared as in Ref. |3j, ro¬ 
tational invariance is realized approximately for bound 
states that fit within the lattice volume, i.e., there is 
an approximate orthogonality between operators sub¬ 
duced from different values of J and M. Using a differ¬ 
ent method to smear operators over many lattice sites, 
Ref. [1], finds that the leading rotational-symmetry vio¬ 
lating contributions from the finite lattice spacing, a, are 
suppressed by a s a 2 as a —> 0, where a s = g 2 /(ATr) and 
g s is the strong coupling constant. 

The analyses of Refs. mm obtain lattice correlation 
matrices that are approximately block-diagonal, with 
each block corresponding to a group of operators sub¬ 
duced from the same value of J. Lattice excited states 
are found to be created predominantly by operators sub¬ 
duced from a single value of J. That identifies the spins 
of excited states created by single-hadron operators. 

However, multiparticle operators must be considered in 
view of the fact that excited states of hadrons show up as 
resonances in the scattering of two hadrons. For example, 
the A baryon, is a wide resonance in nN scattering for 
spin J = | and isospin / = § [5j. The quantum numbers 
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match those of a three-quark state in the quark model, 
or in lattice QCD. 

In a well-known paper, Liischer|6. has shown how 
phase shifts for the scattering of two hadrons can be ex¬ 
tracted from lattice QCD. The Liischer method has been 
extended to systems with nonzero total momentum, un¬ 
equal masses, coupled channels and intrinsic spins la¬ 
ma- It has been applied to coupled-channel scattering in 
Refs. P3HH]- 

In principle, the spins of scattering states can be iden¬ 
tified by matching i), the pattern of approximately de¬ 
generate states found in lattice irreps with ii), the pat¬ 
tern of subductions of a given spin to the irreps of the 
lattice symmetry group. However, the combination of 
accidental degeneracies and statistical uncertainties can 
give ambiguous patterns of lattice states, e.g., the assign¬ 
ment of two or more values of J may be compatible with 
them. 


In this work, sets of lattice operators based on partial 
waves are developed for the scattering of two hadrons 
with zero total momentum, and operators based on he¬ 
licity are developed for nonzero total momentum. The 
operators are designed to facilitate the identification of 
the spin in calculations of phase shifts. 

Periodic partial waves and their orthogonality proper¬ 
ties are discussed in Section [TTJ Operators based on par¬ 
tial waves are defined in Section III Using the partial- 
wave operators, intrinsic spins are included to obtain op¬ 
erators that create lattice states at zero total momentum 
corresponding to total spin J in Sec. m These are called 
parent operators. They are subduced to irreps of Oh for 
use in lattice QCD. The subduced operators are called 
descendant operators. Sets of mutually orthogonal par¬ 
ent operators are identified that correspond to a range of 
different total spins. 

Section IV C develops parent operators based on he¬ 
licity for nonzero total momenta. These are subduced 
to lattice irrep operators, the descendant operators, that 
involve combinations of ±A. m Sets of mutually orthog- 
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onal operators are identified that correspond to a range 
of different helicities. Section[V]provides a brief summary 
of the work. 


II. PERIODIC PARTIAL WAVES 

Periodic boundary conditions in spatial directions ap¬ 
ply to quark propagators on a lattice. Thus, lattice states 
created by quark-field operators generally are periodic 
in spatial directions. In spectroscopy based on single¬ 
hadron operators, the effects of the boundary conditions 
are reduced if the operators create states that fit approx¬ 
imately within the lattice volume. Scattering states are 
affected strongly by periodic boundary conditions. 

Generally, scattering calculations incorporate the peri¬ 
odic boundary conditions by using operators that involve 
combinations of plane-waves with momenta restricted to 
k = , where L is dimension of the lattice in units of 

the lattice spacing and n is a vector with integer compo¬ 
nents along the axes of the lattice. 

Partial waves, /q m (x) = j e (k\x\)Y^ m (x), involving 

momenta k = 27I / | n . n o longer obey the same periodicity 
as the plane waves. However, partial-wave operators cre¬ 
ate lattice states that must obey the periodic boundary 
conditions. Such lattice states can be discussed in terms 
of periodic partial waves that are constructed from the 
formula, 


[/,, m (x)] W =C^/^(x + nL). (1) 

n 

The integer-valued components of n are limited by —N < 
m < N for i = 1 to 3 and the factor, C = (2V+1) -3 , pro¬ 
vides a normalization. In the limit N —> oo, the function 

[,/tm(x)] (per) is periodic, 

[/ f , m (x + nL)] W =[/, m (x)] W , (2) 

where n is any vector with integer components. 

Periodic partial waves for orbital angular momentum, 
£, and projection along the z-direction, m, have restricted 
orthogonality. Consider the overlap matrix, 


O e ,m'-,t,m = J d 3 x [je'(k\x\)Yl m ,(x] (per) 

x[j e (k\x\)Y eim (x)] {per) . (3) 

The integration symbol in Eq. ^ should be read as a 
sum over (L+ l ) 3 lattice points that are inside or on the 
boundaries of the L 3 lattice volume. 

In order to make diagonal overlaps equal to unity, de¬ 
fine normalization factors and construct the normalized 
matrix of overlaps as follows, 




Zv 


1 


,m' :Z,m 


,m' -^£,7 


-O 


i' ,m' :£,rri‘ 


( 4 ) 


TABLE I: Pattern of even-parity matrix elements 
1 for diagonal elements, blank for zero elements and x for large 
off-diagonal elements. The top row shows the l, m indices and 
the left column shows the (!,m! indices. 



0,0 

2,2 

2,1 

2,0 

2,-1 

2,-2 

4,4 

4,3 

4,2 

4,1 

4,0 

4,-1 

4,-2 

4,-3 

4,-4 

0,0 

1 






X 




X 




X 

2,2 


1 




X 



X 




X 



2,1 



1 







X 




X 


2,0 




1 



X 




X 




X 

2,-1 





1 



X 




X 




2,-2 


X 




1 



X 




X 



4,4 

X 



X 



1 




X 




X 

4,3 





X 



1 




X 




4,2 


X 




X 



1 




X 



4,1 



X 







1 




X 


4,0 

X 



X 



X 




1 




X 

4,-1 





X 



X 




1 




4,-2 


X 




X 



X 




1 



4,-3 



X 







X 




1 


4,-4 

X 



X 



X 




X 




1 


Numerical calculations for finite lattice spacings show 
that, 


Zg',rn']£,m — 0 , 

if £’ 7 ^ £ (mod 2) and m! 7 ^ m (mod 4), (5) 

holds to high precision. The pattern of zero elements in 
Eq. § is the same as for Luscher’s matrix, A4e\m'-,e,m, 
that is central to the method for determining phase shifts. 
As shown in Ref. [ 6 ], the zeroes are exact because of 
cubic symmetry for a finite box with no discretization. 
For a discrete lattice spacing, the orthogonality is found 
in Appendix [A] not to depend on the momentum in the 
spherical Bessel functions, or on the difference between 
scattering wave functions and spherical Bessel functions. 
In a lattice QCD calculation, the orthogonality should be 
realized approximately after averaging over gauge config¬ 
urations. 

There is one exception to the mod 2 in Eq. (§: overlaps 
of periodic partial waves involving £’ ,£ = 0 , 2 , or 2 , 0 , are 
zero. Thus, the £ = 0, 1 and 2 periodic partial waves are 
orthogonal to one another, except for the m = ±2 values 
that are affected by the mod 4 rule, i.e., the {£' ,m'\£,m) 
= ( 2 , 2 | 2 , — 2 ) and ( 2 , — 2 | 2 , 2 ) overlaps are of order 1 . 

Periodic partial waves of opposite parities are orthog¬ 
onal. Tables [I] and [id] show the patterns of orthogonality 
of ; £ im for partial waves of the same parity, using 

a blank to indicate a matrix element that is zero, 1 to 
indicate a diagonal element and x to indicate a large off- 
diagonal element. 

Periodic operators with £ = 0, 1 and 2 unavoidably cou¬ 
ple to states of higher £ because of the cubic symmetry. 
For example, a periodic operator for angular momentum 
£ = 1 couples with £ = 3, 5, ■ • • states. Thus, assuming 
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TABLE II: Pattern of odd-parity matrix elements 
Symbols are the same as in Table |T] The top row shows the 
t , m indices and the left column shows the £', m! indices. 



1,1 

1,0 

1,-1 

3,3 

3,2 

3,1 

3,0 

3,-1 

3,-2 

3,-3 

1,1 

1 





X 




X 

1,0 


1 





X 




1,-1 



1 

X 




X 



3,3 



X 

1 




X 



3,2 





1 




X 


3,1 

X 





1 




X 

3,0 


X 





1 




3,-1 



X 

X 




1 



3,-2 





X 




1 


3,-3 

X 





X 




1 


Acting on the vacuum, this operator creates a two-meson 
state with equal and opposite momenta, 

m !(k) 0 )?t 4(—k, 0) 10) = |k,-k)^-. (9) 

In a cubical volume of dimension L, the state is periodic 
when the momentum is k = where n is a vector 

whose components are integers. 

In the sink operator for two scalar mesons, destruction 
operators, a k >, are kept and in general the momentum 
of the sink operator, k', has a different direction but the 
same magnitude as k, 


TOi(k', 0 )m 2 (-k', 0 ) 


dk'CL—k' 


2E k ' 


( 10 ) 


that £ > 5 can be neglected, there is a need to distinguish 
which states are £ = 1 versus i = 3. The blank elements 
in Tables [I] and [IT] show that it is possible to construct 
periodic operators for £ = 3 that are orthogonal to those 
for £ = 1. For example, £,m = 3,2 is orthogonal to £,m 
= 1,1, 1,0 and 1,-1. Thus, including a tracer operator 
that couples to states with £ > 3 but does not couple 
to £ = 1 states, provides a means to distinguish lattice 
states involving £ = 3 from those involving £ = 1. Sim¬ 
ilarly, states involving £ = 4 can be distinguished from 
ones involving £ = 2 by use of suitable tracer operators. 
That provides the motivation for the use of partial-wave 
operators corresponding to distinct values of £ and to. 


Correlation functions involve the operators at t = 0 
and the Euclidean-time propagator, which is written as 
a sum over a complete set of eigenstates, |ra), of the lattice 
Hamiltonian, as follows, 

C(t) = y^(0|mi(k',0)TO 2 (—k', 0)| n)e~ Ent 

n 

x (n\m\ (k, 0 )m 2 (—k, 0 )| 0 ), 

yielding a correlation function for plane-wave states of 
two mesons, 

c (t) = Yl ( k '>- k »e^ nt {n|k,- k ) ( 11 ) 

n k 


III. PLANE-WAVE, LATTICE IRREP AND 
PARTIAL-WAVE OPERATORS 

A. Plane-wave operators 

Consider two point-like fields for distinct, scalar 
mesons that are labeled by indices 1 and 2 . 

The field operator for a single scalar meson at t=0, 
keeping only the creation part, is 

m!(x ' 0) = (JjWF / d3p (^ji7i“l^‘ P ^ (6) 

Fourier transformation of the field gives the creation op¬ 
erator for a plane-wave state, 

™I( k ,0) = ^ 3/2 J d 3 x e ik x m{(x, 0), 

= a ‘ c ( 2 £? k ) 1 / 2 - ^ 

For total momentum and time both equal to 0, a two- 
meson operator is, 

m i (k, 0 )mJ>(—k, 0 ) = 4 at k ^w- ( 8 ) 


B. Lattice irrep operators 

Lattice irrep operators involve sums of plane-wave op¬ 
erators over sets of momenta that are related to one an¬ 
other by rotations that leave the lattice invariant. In 
general, the use of lattice irrep operators leads to mix¬ 
ings of various partial waves. 

An illustrative case is row 2 of the T\ lattice irrep, for 
which the source operator is a sum of two terms, 

7,"4(k,0)r?4(—k,0) - *m{(—k,0)7r4(k,0), (12) 

where k is along the z direction. The sink operator has 
a similar form except with annihilation operators. The 
correlator can be written as, 

C Tl ,2(t) = ^2 E ^ (V’ - k 'l - <~ k '’ k 'l) I") 

k n 

xe _E » t <n|i(|k,-k) — | — k, k)). (13) 

Inserting complete sets of states, |x 1 ,x 2 ), and using x = 
xi — x 2 , X = i(xi +X 2 ) ( similarly for x' and X'), leads 




























to 


Ct ^) = ^2 E e Ent 

k n 

X J <Px' 1 (Px' 2 ^(e i!k '- x ' 

x J d 3 x 1 d 3 x 2 ^ n ( X ,X)^e ik x - e~ ik ^, (14) 


where 'P ra (x, X) = (x x , x 2 1 n) is the wave function of state 
|n). For this example, partial-wave expansion of the fac¬ 
tors involving plane-waves contain odd i values, 

Ee lk x - e“ lk x ) = ^i e je(k\x\)Yl 0 (k) Y e , 0 (x). 

odd e 

(15) 

Thus, the correlator becomes, 

C Tl ,2(t) =E e_E ” t i( E 47r/^ >0 (fc , )^^,o(fc)) 

n k odd V 

x(E (*)*2o(*))> (!6) 

odd t 

where 

= J d 3 x j<(A;|x|)y^ m (i)$ n (x), (17) 

and 

$ n (x) = J d 3 X* n ( X ,X). ( 18 ) 

The matrix elements, 4> n ^ jm (fc), show that the coupling 
of operators to lattice states involves the j^{k\x\)Y^ m {x) 
part of the operator. 

When the mesons do not interact, the wave function, 
4> n (x), corresponds to two mesons with the same mo¬ 
menta as in the operator that creates the state. The wave 
function can be expanded in partial waves as follows, 


$n(x) = E ac ’ M ( 19 ) 

C,M 

Substituting this expansion into Eq. ( p~7| ) and taking into 
account the periodic boundary conditions leads to sums 
over the the overlaps of Eq. @. 

When the two mesons interact, the scattering wave 
function generally differs from jc(k |x|). Appendix [A] de¬ 
scribes numerical tests that show that the interactions 
do not have a significant effect on the orthogonality. 


This is a form that is useful for lattice QCD. In order 
to restrict the source or sink operator to a single partial 
wave, the plane-wave factor is replaced by a partial-wave 
factor, defined as follows, 


_ik-x 


= je(k\-x\)Ye tm (x), 


( 21 ) 


where a factor 47r^Y^ m (fc) is omitted. The partial-wave 
operator is based on the factors that occur in the ma¬ 
trix elements between operators and lattice states, as 

uivalent to projecting the 


in Eq. |17|). The result is ecn 

plane-wave factor e k x in Eq. (20|) to a partial wave using 

Y e>m (k). 

The partial-wave operator is defined by, 


m^(k, 0)m^(—k, 0) 

J 

^ d 3 x i f d 3 x 2 


(2t r) 3 


„zk-: 


J L‘ 


TO' l '(xi)m t (x 2 ). (22) 


When used in a lattice QCD calculation, a partial-wave 
operator creates periodic states. 

The quantization of momentum for noninteracting 
plane-waves that obey periodic boundary conditions im¬ 
plies that q = = |n|, where n is a vector with in¬ 

teger components. Interactions modify the wave func¬ 
tions but they should remain close to the noninteracting 
waves with momenta corresponding to n 2 = 1, 2, 3, 4 • • •. 
While there is no restriction that limits the momentum 
values, the noninteracting momenta provide a reasonable 
choice for the momenta in the partial-wave operators. 
The jc(k |x|) factor can also be viewed as simply provid¬ 
ing a spatial smearing with a smearing parameter, fc, that 
may be chosen to maximize the coupling of the operator 
to states of interest. 

A matrix of correlation functions could be obtained 
by using partial-wave operators with orbital angular mo¬ 
mentum indices (! and m! at the sink point together with 
indices £ and m at the source point, which gives, 


= Y e Ent ^n,t',rn'{k)^ ntm {k), (23) 


where & n ,e,m{k) is defined by Eq. {171. In this case, one 
can select operators corresponding to an orthogonal set of 
states, including tracer operators to indicate the presence 
of higher spins, in order to facilitate the identification of 
spins. However, the partial-wave operators are not used 
directly. They serve as parent operators that are sub- 
duced to lattice irreps before use in lattice calculations. 


C. Lattice partial-wave operators 


In the analysis above, the two-meson operator at t = 0 
contains the plane-wave factor, e* k x , and can be written 
as, 


m^(k, 0)m^(—k, 0) = 

——^ J d 3 x i J d 3 x 2 e tk x to^(xi, 0 )m^(x 2 ,0). (20) 


IV. OPERATORS FOR TWO-HADRONS WITH 
INTRINSIC SPINS 

In lattice QCD, the representations of the intrinsic 
spins of hadrons are orthogonal, like the spin states in 
quantum mechanics, 
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where s is the spin and m s is its projection along the z- 
azis. Parent operators for the scattering of two hadrons 
with zero total momentum are constructed in Sec. IIV Al 
by coupling the intrinsic spin angular momenta with a 
partial wave of orbital angular momentum to make to¬ 
tal spin J and projection M. The final step, outlined in 
Sec. |IVB| is to produce descendant operators by subduc- 
ing the parent operators from J and M to lattice irreps 
of O h . 

When total momentum is nonzero, operators are con¬ 
structed based on total helicity in Sec. IV C[ Lattice irrep 


operators are constructed in Sec. |IVE 


A. Parent operators for total momentum P = 0 


Consider zero total momentum: P = ki + k 2 = 0. 
The relative momentum is k = |(ki — k 2 ) and it follows 
that ki = k and k 2 = —k. The masses of the hadrons 
are mi and m 2 , and the energies are Fj = i/k 2 + m\ 
and E 2 = -\/k 2 + m|. For fixed energies of the hadrons, 
relative momentum, k, is a vector of fixed magnitude, k. 

Two equivalent ways of writing the lattice coordinates 
of the operators are, 

xi = X + -x, x 2 = X- -x, (25) 

and 

X=i(xi+x 2 ), x = Xi-x 2 . (26) 


We switch between them without further comment. 

Consider a product of hadron operator, A Slirris ,(xi,t), 
and hadron operator, B S2t7Usr> (X 2 ,t), where Si and s 2 are 
intrinsic-spin labels, and m Sl and m S2 are labels for the 
z-components of intrinsic spins. The desired parent oper¬ 
ator involves coupling intrinsic spins and relative orbital 
angular momentum to total angular momentum, J, and 
projection, M, as follows, 


(AB)jp M (P = 0,k,t) =f d 3 x±f d 3 x 2 C 

J ” «Ji ftn 

A Slt m B1 ( x l) t)B S2 ,m a2 ( x 2, t)> 


JM 


Sl,«2 
m sl , m S2 


(27) 


where Cjm denotes the following product of the usual 
Clebsch-Gordan coefficients for coupling the spins and 
orbital angular momenta, 

Cjm = (s 1 m Sl s 2 m S2 \Sms}(Sms(im\ JM), (28) 


and the partial wave of relative momentum is given by 

Eq. @. 

A matrix of correlation functions for the hadron- 
hadron system based on parent operators is, 


Cj'm'-.jm (P = 0, M) = <0| ( AB)j, m , (P = 0, k,t) 
x JAB) jm (P = 0, k, 0) |0>. (29) 


TABLE III: Mutually-orthogonal parent operators for the 
scattering of two hadrons with total intrinsic spin 5 = 0, for 
example, 7T7 t or NN. The final column provides conventional 
labels of the form 2S+1 Lj for the NN system with 5 = 0. 


Operator 

J P ,M 

i,m 

25+1 L j 

1 

0 + ,0 

0,0 

1 So 

2 


1,1 

x Pi 

3 

f,o 

1,0 

'Pi 

4 

2+,2 

2,2 

1 d 2 

5 

2+1 

2,1 

'd 2 

6 

2+,0 

2,0 

1 d 2 

7 

3",3 

3,3 

] f 3 

8 

3",2 

3,2 

'Fa 

9 

4",3 

4,3 

x g 4 


Sets of mutually-orthogonal parent operators for two- 
hadron scattering states ar e giv en in Table [TTT| for total in¬ 
trinsic spin S = 0, in Table IV for S = \ and Table |v| for 
5=1. The tables provide the quantum numbers corre¬ 
sponding to nonzero Clebsch-Gordan coefficients for the 
coupling of intrinsic spins with orbital angular momenta 
t = 0, 1, 2, 3 and 4. Some of the listed operators are or¬ 
thogonal to others because their couplings involve pairs 
of £, m values that are orthogonal, or pairs of m s values 
that are orthogonal. Other operators, for example § + | 
and | |, involve combinations of Clebsch-Gordan coeffi¬ 
cients that give orthogonality. Note that each £, m com¬ 
ponent involved must be normalized to the same value 
in order to realize the orthogonality due to the Clebsch- 
Gordan coefficients. I 1 ! 


It is straightforward to add operators for other mo¬ 
menta in order to build a variational basis of operators 
with respect to energy within each set of J, M values. 
For each additional momentum included, all operators of 
Table [TTT] ( or Table IV or Table |v|) can be used. Orthog¬ 
onality is then between subsets of operators that have 
different J, M values. 


Isospin and the Pauli Principle provide limits on the 
partial wave operators that should be used. For exam¬ 
ple, the 7T7T system can occur in states of isospin 0, 1 or 
2 and the system must be symmetric with regard to the 
exchange of the pions. Similarly, the NN system can 
occur in states of isospin 0 or 1 and the system must be 
antisymmetric with regard to exchange of the two nucle¬ 
ons. Operators that do not obey the correct symmetry 
should be omitted. 


[1] It is straightforward to use the orthogonality properties of C-G 
coefficients to project out the i, m components corresponding to 
a given m s . 
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TABLE IV: Mutually-orthogonal parent operators for the 
scattering of two hadrons with total intrinsic spin 5 = for 
example, nN. Each row of the table gives the combinations 
of £, m and m s quantum numbers that occur in the Clebsch- 
Gordan expansion to form J p and M. 


Operator 

J p M 

£, m m s 

£, m m s 

1 

1 + 1 

2 2 

0,0 i 


2 

1 - 1 

2 2 

1,0 i 

1,1 -i 

3 

1 - 1 

2 "2 

1,-1 i 

1,0 -i 

4 

3 - 3 

2 2 

1,1 i 


5 

3 - 1 

2 2 

1,0 i 

1,1 -i 

6 

3 + 3 

2 2 

2,1 i 

2,2 -I 

7 

3 + 1 

2 2 

2,0 i 

2,1 -i 

8 

5 + 5 

2 2 

2,2 i 


9 

5 + 3 

2 2 

2,1 | 

2,2 -i 

10 

5 + 1 

2 2 

2,0 i 

2,1 -I 

11 

5 - 5 

2 2 

3,2 | 

3,3 -i 

12 

7 “ 5 

2 2 

3,2 i 

3,3 

13 

7 + 7 

2 2 

4,3 i 

4,4 


TABLE V: Mutually-orthogonal parent operators for the 
scattering of two hadrons with total intrinsic spin 5=1, 
for example, NN. Each row of the table gives the combi¬ 
nations of £, m and m s quantum numbers that have nonzero 
coefficients in the Clebsch-Gordan expansion to form J and 
M. The final column provides conventional labels of the form 
2 S+ 1 Lj for the NN system with 5=1. 


Operator 

J 

M 

l, m 

m s 

l, m 

m s 

£, m 

m s 

2s+1 Lj 

1 

0 

1 

0,0 

i 



3 Si 

2 

0 

0 

0,0 

0 



3 Si 

3 

0 

-1 

0,0 

-1 



3 Si 

4 

0 

0 

1,1 

-1 

1,0 

0 

1,-1 

i 

3 Po 

5 

1 

1 

1,1 

0 

1,0 

1 


3 Pi 

6 

1 

0 

1,1 

-1 

1,-1 

1 


3 Pi 

7 

2 

2 

1,1 

1 



3 P 2 

8 

2 

1 

1,1 

0 

1,0 

1 


3 p 2 

9 

2 

0 

1,1 

-1 

1,0 

0 

1,-1 

i 

3 p 2 

10 

1 

1 

2,2 

-1 

2,1 

0 

2,0 

i 

3 Di 

11 

1 

0 

2,1 

-1 

2,-1 

1 


3 Di 

12 

2 

2 

2,2 

0 

2,1 

1 


3 d 2 

13 

2 

1 

2,2 

-1 

2,1 

0 

2,0 

i 

3 d 2 

14 

2 

0 

2,1 

-1 

2,-1 

1 


3 d 2 

15 

3 

3 

2,2 

1 



3 d 3 

16 

3 

2 

2,2 

0 

2,1 

1 


3 d 3 

17 

3 

1 

2,2 

-1 

2,0 

1 

2,0 

i 

3 d 3 


B. Descendant operators in lattice irreps for P = 0 


When the total momentum of a two-hadron operator is 
zero, the lattice is invariant with respect to the elements 
of the octahedral group. The parent operators given 


above are subduced to irreps of the octahedral group, 
which produces descendant operators that can be used 
in lattice calculations. 


The subduction coefficients, S '^^, are given in Ref. [T] 
for integer spins and in Ref. [2] for half-integer spins. 
They should be applied as follows to the operators of 
Eq. (271 that are labeled by J,M to obtain operators 
that are labeled by the octahedral group irrep A and row 


{AB) Ar (ki, k 2 , t) = J2 S a™ ( ab )jm (ki, k2, *)■ (30) 

M 


As shown in appendix B of Ref. [2], the subduction 
coefficients are real and they obey, 


6 A',r' 6 A,r = <W,A<V,r- 


(31) 


M 


This property causes descendant states to be orthogonal 
with respect to J when they are obtained by subducing a 
set of mutually orthogonal parent states. This is shown 
for quantum states of angular momentum, | J, M ), and 
their descendant states, labelled by lattice irreps, that 
are subduced from a single value of J, i.e., |A[J],r) = 
E m s a I J , M )» as follows, 


(A'[J'],r'|A[J],r) = £ S0'(J 1 ,M'\J,M), 


M,M' 


E qJ',M' c J,M c x 

<>A',r' ‘-’A ,r d J',J d M',M 


M,M' 


= $A',A Sr',rSj',J- 


(32) 


Spin identification is based on an orthogonal set of 
parent operators. A matrix of correlation functions is 
calculated based on the descendant, irreducible opera¬ 
tors. Because of the orthogonality of Eq. (32), the ma¬ 
trices are approximately block-diagonal with respect to 
the values of the parent spin, J , from which the irre¬ 
ducible operators are subduced. Use of the variational 
methodjT8l T9] to extract the energies of states, and op¬ 
erator overlaps M, Znm = (n|O m |0), to identify the 
couplings of operators to the states, then identifies which 
parent spins, J, are predominant in the creation of a lat¬ 
tice state. Generally, when a spin-J operator couples to 
a state, that identifies J as a lower bound of the spin. 
If two or more operators couple to a state, the highest 
value of J identifies the lower bound of the spin. 


The phase shift can be calculated using the Luscher 
method, which is based on the deviation of the energy 
of a scattering state from the energy when there is zero 
interaction. The latter energy may be calculated from 
the masses of the hadrons and the quantized momenta 
that obey periodic boundary conditions. 
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C. Two-hadron operators with total momentum 

P A 0 


The matrix of correlation functions for the two-hadron 
system based on parent operators is, 


Nonzero total momentum is of considerable interest 
because it provides a way to increase the number of c.m. 
energies that can be obtained using a given lattice. 

The helicity, A, is the eigenvalue of J z in the c.m. 
frame. Invariance with respect to boosts along the di¬ 
rection of the total momentum shows that helicity takes 
the same value in the c.m. frame and the lattice frame, 
where it can be determined. 

Because J > |A|, the use of operators for a range of 
helicity values provides the following signature. A state 
with spin J couples to the operators with A < J, and not 
to the operators for higher helicity. 

It is useful to decompose coordinates and momenta 
into parts that are perpendicular and parallel to P, which 
are denoted by subscripts, T and ||, respectively. Rota¬ 
tions that keep P invariant also leave the parallel com¬ 
ponents of other vectors invariant, e.g., Ami, k 2 \\ and 

fcy = — /c 2 1 ^. Vectors that are perpendicular to 

the total momentum, e.g., ku_ = — k 2 x = kx, are in¬ 
variant with respect to boosts along the direction of the 



and E 2 = ^jm\ + k\ + k^, where k± = |kjJ is invari¬ 
ant with respect to rotations that keep P invariant. 

D. Parent helicity operators for P A 0 


C\'-x{P,k\i,k ± ,t) = (0| (AB) a , (P, k\\,k x ,t) 
x(AB) A (P,fc||,fc x ,0)|0). (36) 

E. Descendant operators in lattice irreps for P A 0 

When the total momentum of the two-hadron system 
is nonzero, the lattice symmetry group is C nv (alterna¬ 
tively known as Dic„ for the double group), for n =4, 3, 
or 2. For C 4 ,,. the total momentum is directed along one 
of the axes of the cube; for C 3v , it is along a body diag¬ 
onal; for C 2v it is parallel to a face diagonal. For each 
of these cases, we label the direction of total momentum 
as the z-axis. Altmann and Herzig [20] have developed 
the subductions from SU (2) to irreps of the point groups 
C nv \ their notation is followed in this work. Equivalent 
results are given in Refs. D2UZHE2I. 

The C nv group applies for integer spins and the C nv 
double group applies for half-integer spins. The group 
elements are n rotations about the z direction, R z (9 ), 
with equally spaced 9 values in the range 0 < 9 < 2n 
(or 2 n rotations in the range 0 < 9 < 47r for the double 
group). In addition, there is an equal number of reflec¬ 
tions, R z (9)XV yz , that involve the same rotations and a 
reference reflection, which is chosen to be II yz , the re¬ 
flection in the yz- plane. The group elements transform 
helicity states as follows, 


Two-hadron operators are labelled by the total helic¬ 
ity when P A 0. For two spin-zero hadrons, the helicity 
values are constructed in the lattice frame using a cylin¬ 
drical function in place of the spherical harmonic that 


appears in Eq. (211, as follows 


0 ikx ‘X ± 


= J m (k±x ± )e im ** 


(33) 


where (j) x is the azimuthal angle of xj_ and J m is a Bessel 
function. An overall factor j s omitted, where 

(f>k is the azimuthal angle of kx, because it plays no role 
in the final analysis. 

A two-hadron parent operator incorporates intrinsic 
spins of the hadrons as follows, 


(AB) a (P, k\\ ,k±,t) = f d 3 xi j d 3 x 2 e* p ' x e ifc n x ii ^ 

J ” fit . 


Sl,S2 

m si ,m S2 ,tti 


C\ 


_,ik_L -X_L 


A si ,m si 


(xi ,t)B S 2 l m.Ax 2 ,t), (34) 


where C\ denotes the following product of a Clebsch- 
Gordan coefficient for spin and a delta function, 


C\ = ( s 1 m Sl s 2 m S 2 \Sms) S mg+m, A* (35) 

This operator is constructed to create a state with helic¬ 
ity, A, that is the sum of the orbital helicity, m, and spin 
helicity, ms- 


R Z (9)\X) = e~ lXe \X), 

n„*|A)=»?|-A), (37) 

where rj = (—1 ) J P and P is the parity. The phase factor, 
rj, is equal to f] in the analysis of Ref. SU, where it is 
shown to follow from the choice of II yz as the reference 
reflection, 

Representations of C nv involve a linear combination 
of helicity states, |A) and | — A). In this work, they 
are chosen to diagonalize the reference reflection and are 
labeled by the positive helicity, A > 0, and a sign, ±, 
which serves as the row index of the representation, as 
follows, 


|^>±) = V^( l A ) ±r ?l _A ))- ( 38 ) 

The action of group elements is to transform the |A, ±) 
representations into linear combinations of one another, 

R z (d)|A,+) = cos(A0)| A, +) — A, —), 

R z (9) | A, —) = — isin(X9) | A, +) + cos(X9) | A, — ), 
n yz |A,±) = ±(—1)' 7 | A, ±), (39) 

where the eigenvalue of the reference reflection is seen to 
be ±(—1) J . Note that a two-dimensional representation 










is reducible to two one-dimensional irreps when sin(Xd) 
vanishes for each group rotation angle, 0. Otherwise 
|A,+) and | A, —> provide the rows of a two-dimensional 
irrep. 

Descendant operators for P / 0 transform as irreps of 
C nv . They involve the positive helicity, A > 0, and the 
row index, ±, and have a composition similar to Eq. (38), 


(AB)x,±(P,k\\,k±,t) 



P,k\\,k±,t) 


±r 1 (AB)_ x (P,k ll ,k ± ,t)). (40) 


The rows of a two-dimensional irrep are orthogonal 
and they produce the same correlation function as the 
average over rows, which is the same as the average over 
± signs. Using the descendant operators of Eq. (40), the 
average over rows is, 

= \ (Ca,+;A,+ + CV;A,-) • (41) 

A little algebra shows that this is equal to the aver¬ 
age over the helicities, A and —A, that appear when the 
expression above is expanded in terms of the correlation 


functions Ca',a of Eq. (36). Using 77*77 = 1, the row av¬ 
erage is evaluated as follows, 


C A 


— ^ (Ua,A + T]*C— A,A + + 77* 77C — A,—a) 

+ ^ (Ca,a, - V*C-x,x - r]Cx,-x + rfr)C~x,-x ), 


= 2 (C x ,x + C-x-x) • 


(42) 


The parity-dependent phase factor, 77 , cancels out of the 
average over rows. 

Each one-dimensional irrep also involves both A and 


—A, as in Eq. (38). For most purposes, it is adequate 


to average the two one-dimensional irreps for a given 
helicity, A, in the same manner as the rows of a two- 
dimensional irrep. When this is done, the result is the 
same as in Eq. (421, i.e., the parity-dependent phase fac¬ 


tor, 77 , cancels out. Thus, correlation functions can be 
calculated directly from parent operators, as in Eq. (36), 


and then averaged over the ±A values in order to obtain 
the same result as would be obtained from averaging over 
the rows of correlation functions based on the irrep op¬ 
erators of Eq. diOl. 


Although the average over helicities provides a sim¬ 
ple and straightforward analysis, the irrep operators are 
given in Appendix [B] for completeness. 

The orthogonality of periodic orbital-helicity functions 
involves the following matrix of overlaps, 

Zm'm = J d 3 x[j m (k ± X ± )e im **Y Per) ' X 

[Jm'(k±x ± )e im '^] {per) . (43) 


TABLE VI: Pattern of matrix elements Z m / rn for group Ci v . 
For group Cz v , the x’s shift to horizontal positions that are a 
multiple of three units away from the diagonal and for group 
Civ the x’s move to positions that are a multiple of 2 units 
away from the diagonal. Symbols have the same meaning 
as in Table |I] The top row shows the m indices and the left 
column shows the nn! indices. 



-4 

-3 

-2 

-1 

0 

1 

2 

3 

4 


-4 

1 




X 




X 


-3 


1 




X 





-2 



1 




X 




-1 




1 




X 



0 

X 




1 




X 


1 


X 




1 





2 



X 




1 




3 




X 




1 



4 

X 




X 




1 



TABLE VII: Mutually orthogonal helicity operators for group 
Ci v and total intrinsic spin S = 0. Each row lists the irrep 
label, A, and the m values included in the operator. For 
groups Czv and C-iv ., only rows 1 and 2 involving m = 0 and 
±1 are orthogonal. 


No. 

A 

m 

-m 

1 

0 

0 

0 

2 

1 

1 

-1 

3 

2 

2 

-2 


This matrix is normalized as in Eq. Q in order to make 
diagonal elements equal to one. Then the zero elements 
depend on index n of group C nv as follows, 


Z m ’ ,m =0 if m' 7 ^ 771 ( mod n). (44) 


Table JVT] shows the pattern of matrix elements based on 

Eq. @7 

The best case for spin identification is group C± v be¬ 
cause it has the largest range of orbital helicities that 
can be included in an orthogonal set. The orbital-helicity 
states that are orthogonal are limited to 777 = 0 and ±1 for 
both C% v and C 2 V ■ For each of these groups, the result¬ 
ing sets of mutually-orthogonal operators for two-hadron 
scattering states are g iven i n Table [VTT| for total i ntrin sic 


spin S = 0, in Table 


VIII 


for S' = ^ and Table 


IX 


for 


S = 1. As each irrep operator includes helicity values, 
A and —A, quantum numbers m and —m are included 
for the cylindrical harmonics. Orthogonality for S = 0 is 
limited because each 777 -component must be orthogonal 
to all the other 777-components in the table. For example, 
?77 = ±3 cannot be included because 777 = ±1 is included. 


When intrinsic spin states are included, orthogonality 
depends on each m, m s pair being orthogonal to each 
other pair. The orthogonality of spins allows more oper¬ 
ators to be included. 
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TABLE VIII: Mutually orthogonal helicity operators for 
group C 4v and total intrinsic spin S = Each row lists 
the irrep label, A, and the m, m s values included in the oper¬ 
ator. For groups C 3 V and C 2 v, only the rows involving m = 
0 and ±1 are orthogonal. 


No. 

A 

m ,m s 

-m,-ra s 

1 

1 

2 


0,-1 

2 

1 

2 


-1,1 

3 

3 

2 

1 ,1 

- 1,4 

4 

5 

2 

2,1 

- 2,4 


TABLE IX: Mutually orthogonal helicity operators for group 
Ci v and total intrinsic S = I. Each row lists the irrep label, 
A, and the m, m a values included in the operator. For groups 
Czv and C 2v , only the rows involving m = 0 and ±1 are 
orthogonal. 


No. 

A 

m, m s 

-m, —m s 

1 

0 

0,0 

0,0 

2 

0 

1,-1 

-1,1 

3 

1 

0,1 

0,-1 

4 

1 

1,0 

-1,0 

5 

1 

2,-1 

-2,1 

6 

2 

1,1 

-1,-1 

7 

2 

2,0 

-2,0 


V. SUMMARY 

Operators are developed to facilitate the identification 
of spins in lattice QCD calculations of the scattering of 
two hadrons. In the case that the total momentum of 
the hadrons is P = 0, partial-wave operators are used to 
create periodic partial waves of orbital angular momen¬ 
tum, £, and projection, to. The operators are based on 
momenta of magnitude k = —j where is n is a vec¬ 
tor with integer components. These are the momenta 
for noninteracting plane waves on the cubic lattice, how¬ 
ever, other values of k could be used to provide a differ¬ 
ent radial smearing of the operators. Phase shifts may 
be obtained using the Liischer method and its extensions. 
The periodic partial waves are orthogonal with respect to 
i and m, subject to restrictions imposed by the cubical 
lattice. The partial waves provide an advantage relative 
to the use of combinations of plane waves in that there is 
freedom to select sets of orthogonal operators based on 
specific £ and to values. 

The essential point of this work is to differentiate be¬ 
tween lattice irreps that are subduced from different con¬ 
tinuum spins. For example, the use of mutually orthog¬ 
onal £ = 1 and £ = 3 parent operators based on par¬ 
tial waves provides orthogonal descendant operators af¬ 
ter subduction to the Xj lattice irrep. Strong coupling of 
a lattice state to one of these operators but not the other 
provides a signal for the spin of the state. While contri¬ 
butions of £ > 5 can occur in the subduced operators, the 


short range of hadronic interactions and the centrifugal 
barrier at high £ suggest that the role of higher spins is 
suppressed. 

Parent operators for total spin, J < 4, and projec¬ 
tion, M, are obtained by coupling the intrinsic spins of 
the hadrons and the partial waves of orbital angular mo¬ 
mentum. Sets of mutually-orthogonal parent operators 
corresponding to different values of J and M are given 
for total intrinsic spin S = 0,1 and 1 based on omitting 
operators that are not orthogonal because of the cubic 
symmetry. When these parent operators are subduced to 
lattice irreps, their descendants retain orthogonality and 
lead to matrices of correlation functions that are block 
diagonal with respect to the different values of J. That 
allows the spin of a scattering state to be identified as the 
spin of the parent operators that dominate the creation 
of the state. 

In the case that the total momentum is P ^ 0, parent 
operators are designed to create states of total helicity 
by combining spin helicities and orbital helicity based 
on cylindrical harmonics. Mutually-orthogonal sets of 
parent operators for a sufficient range of total helicity 
values can be used to facilitate the identification of the 
spin, J, of a state. The signature is that a state of spin 
J couples to operators for all total helicity values A < J, 
and not to operators for A > J. However, there is a 
limitation owing to cubic symmetry. Total momenta that 
are parallel to an axis of the cubical lattice (symmetry 
group Ci v ) allow the widest range of helicity values and, 
thus, provide the best case for spin identification based 
on helicity. 

For either the partial-wave operators or the helicity 
operators, a variational basis can be constructed by us¬ 
ing additional copies of the sets of mutually-orthogonal 
operators based on different momentum values. 
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Appendix A: Lattice orthogonality of periodic 
partial waves. 

A generalization of Eqs. © and Q is considered in 
which the lattice state is assumed to be a scattering state 
described by the Liischer formalism for various values of 
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the phase shift, 


Zt'Mm = N 1 N 

w Z‘,m'M i,m 

X f d 3 x [j t (fc |x| )yj m , (£’] {per) (|x| )^,m (®)] {per) , 

(Al) 

where \E^(|x|) is the radial wave function of the scattering 
state. When the phase shift, 5g is zero, ^(|x|) reduces 
to a spherical Bessel function, j/(fc|x|). 

According to the Liischer formalism, a scattering state 
for angular momentum £ and phase shift Se involves a 
momentum value, k ext , in the exterior region, r > R, 
where R is the range of interaction. Momentum k ex t 
is relat ed to the e ner gy of the scattering state through 
W = y/rnl + k1 xt + \Jm\ + k^ xt , as shown by Liischer |6]. 
The value of k ext also is related to the phase shift as 
follows, 


6 e= jn - (j)(q), (A2) 

where the function <j>{q) is obtained from 


cot{4>) 


Zoo(l;q 2 ) _ 1 y- 1 ani 

n 3 / 2 q 2ir 2 q ^ (n 2 - q 2 )' 1 1 


The integer j is chosen so that <p(q) is a continuous func¬ 
tion of q and obeys the boundary condition <p(0) = 0 . 
These formulas neglect contributions from £ = 4 or higher 
partial waves. 

Since the interaction is not known, it is necessary to 
use a model to approximate the scattering wave function 
in the interior region, 0 < r < R. A square-well potential 
gives a momentum, kt n t, in the interior region and the 
scattering wave function is, 


^e{r) = jiihntr), 0 <r < R, 

T/>(r) = A[cosSeje{k ext r) - sin8(.nt{k ext r)\,R < r < L/ 2 , 

(A4) 


where jt and rig are spherical Bessel functions. 

The scattering wave function is then obtained by 
choosing a phase shift, finding k ex t from Eq. (A2| and 
matching the interior and exterior wave functions in 
Eq. (A4), and their derivatives, at r = R to obtain mo¬ 
mentum ki n t and coefficient A. 

Calculations have used R = L/4 as the range of the 
interaction and Si/n = [0.0,0.1,0.25,0.5,0.75,0.9]. The 
lowest values of k ex t are used. The cubical volumes con¬ 
sidered are 16 3 and 32 3 . The construction of periodic 
partial waves was done using N = 5, 7, 9 and 11 for the 
sum in Eq.|l]). For all cases, the results confirm that ma¬ 
trix elements shown by blank entries in Tables [I] and [II 


have numerical values less than 10“ 15 in magnitude for 
each considered value of the phase shift. 

Appendix B: Irreducible representations for helicity 
states 


As noted in the text, the two-dimensional representa¬ 
tion of Eq. (38) is reducible when sin(Xd) vanishes for all 
group rotations. Then, the | A, ±) states are not mixed 
by the group rotations, as may be seen from Eq. (39). 


The reduction occurs when X = 0, n/2, n, • • •, where 
n is the index of group C nv . Following the Mulliken 
convention[251 for labelling irreps, the helicity states 
| A, ±) provide irreps as listed in Table |x| For A = 0 or n, 
the | A, +) state transforms as the A\ irrep and the | A, —) 
state transforms as the A 2 irrep. For A = n/2 or 3n/2, 
|A,+) transforms as the B 1 irrep and | A, —) transforms 
as the B 2 irrep. 


For A = 1, the | A, ±) are mixed by the action of group 
rotations for groups C± v and C?, v , but not for C^. For A 
= 3, the rotations mix the | A, ±) for group C± v , providing 
two rows of the E 2 two-dimensional irrep. For groups C^ v 
and C 2 V , there is no mixig and the | A, ±) provide one- 
dimensional irreps A\ and A 2 for C% v and B\ and Bo for 

c 2v . 


For half-integer spins and helicities, the irreps are two- 
dimensional and are given labels Ei , £3 and E, as shown 
in Table |x| An equivalent table may be found in Ref. [22] , 


TABLE X: Helicity irreps for C nv groups, n = ±(—1) J is the 
eigenvalue of the reference reflection and n/(— 1) J is the row 
index, namely, ±. 


A 

n/(-i) J 

Cav 

c 3v 

c 2v 

0 

+ 

Ai 

Ai 

Ai 

0 

- 

A2 

A2 

A2 

1 

2 

± 

Ei 

2 

Ei 

2 

E 

1 

+ 

E2 

E2 

Bi 

1 

- 

Eo 

E ‘2 

b 2 

3 

2 

± 

E3 

B\ ® B 2 

E 

2 

+ 

Si 

e 2 

Ai 

2 

- 

b 2 

E2 

A2 

5 

2 

± 

Ei 

Ai © Bi 

E 

3 

+ 

E2 

Ai 

Bi 

3 

- 

Eo 

A2 

b 2 

7 

2 

± 

Ei 

2 

Ei 

2 

E 

4 

+ 

Ai 

e 2 

Ai 

4 

- 

A -2 

e 2 

A2 
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